Exact results and new insights for models denned over small world networks. 
First and second order phase transitions. II: Applications 
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We apply a novel method (presented in part I) to solve several small- world models for which the 
method can be applied analytically: the Viana-Bray model (which can be seen as a or infinite 
dimensional small- world model) , the one-dimensional chain small- world model, and the small- world 
spherical model in generic dimension. In particular, we analyze in detail the one-dimensional chain 
small-world model with negative short-range coupling showing that in this case, besides a second- 
order spin glass phase transition, there are two critical temperatures corresponding to first- or 
second-order phase transitions. 
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I. INTRODUCTION 

This paper represents the continuation of another pa- 
per [l[ - in the following referred as part I - in which a 
general method to analyze small-world models has been 
presented. Here (part II) we apply it to some cases of 
interest which can be fully faced analytically. Even if the 
exactness of this method is limited to the paramagnetic 
regions, it gives the exact critical behavior and the exact 
critical surfaces, and provide a clear and immediate (also 
in terms of calculation) insight of the physics. As stressed 
in part I, the underlying structure of the non random part 
of the model, i.e., the set of spins staying in a given lattice 
Cq of dimension do and interacting through a fixed cou- 
pling J , is exactly taken into account. When J > 0, the 
small- world effect gives rise to the known fact that a sec- 
ond order phase transition takes place, independently of 
the dimension do and of the added random connectivity 
c. However, when Jo < 0, a completely different scenario 
emerges and - besides a second-order spin glass transition 
- for a sufficiently large c, multiple first and second-order 
ferromagnetic phase transitions may take place. Here we 
will emphasize above all this aspect which, to the best of 
our knowledge, until now was observed only in |2| and in 
the context of small- world neural networks in [31], and in 
our opinion represents a new paradigm in view of inter- 
esting applications in real small-world networks. 

The main advantage of our method lies in its great sim- 
plicity. In fact, what is required to apply it, is to solve 
not the small- world model (a random model) defined over 
Co, but a corresponding non random model still defined 
over Co and immersed in an external magnetic field. This 
implies that we are able to solve analytically small-world 
models whose underlying lattice Co has dimension 0, 1 or 
infinite, as for an ensemble of non interacting units, the 
one dimensional chain, and the spherical model, respec- 
tively. In fact, in all these cases the non random model 
can be exactly solved even when immersed in an external 
magnetic field. 

The paper is organized as follows. In Sees. II and 



III we recall the definition of the small-world models and 
the method, which mainly consists in finding the solution 
of the self-consistent equation () 1 2(1 minimizing the effec- 
tive free energy (|24|) . In Sec. IV we analyze two simple 
cases in which Co has dimension do — 0: an ensemble of 
non interacting units. The small world-model defined on 
it gives rise in particular to the well known Viana-Bray 
model (sometime called random Bethe lattice model) . In 
Sec. V we consider the case in which Co is the one dimen- 
sional chain. Sec. VI is devoted to the spherical model 
- an infinite-dimensional model - defined for arbitrary do 
Finally, in Sec. VII some conclusions are drawn. 



II. SMALL WORLD MODELS 

We consider random Ising models constructed by 
super-imposing random graphs with finite average con- 
nectivity onto some given lattice Co whose set of bonds 
and dimension will be indicated by To and do, re- 
spectively. Given an Ising model - shortly the unperturbed 
model - of N spins coupled over Co through a coupling 
Jo and with Hamiltonian 



Ho — —Ji 



2J ViO-j - h a i ) 



(1) 



and given an ensemble C of unconstrained random graphs 
c, c G C, whose bonds are determined by the adjacency 
matrix elements c,-, 7 - = 0,1, we define the correspond- 
ing small-world model - shortly the random model - as 
described by the following Hamiltonian 



H c ,j = H - V c 



(2) 



the free energy F and the averages (O) 1 being defined in 
the usual (quenched) way as 



(3F 



c£C 



P(c) / dV({J hJ }) log (Z c ,j), 



(3) 
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and 



III. AN EFFECTIVE FIELD THEORY 



(OY = P(c) [ dV ({Ji,j}) (0)\ 1 = 1,2 (4) 

cGC J 

where Z c -j is the partition function of the quenched 
system 



(5) 



(0) C ;j the Boltzmann-average of the quenched system 
(note that (0) c -j depends on the given realization of the 
J's and of c: (O) = {0)c-,J\ f° r shortness we will often 
omit to write these dependencies) 



-PS B;J ({<Ti}) 



Z r ,. 



(6) 



and dV {{Ji,j}) and P(c) are two product measures 
given in terms of two normalized measures dfi(Ji^) > 
and p(cij) > 0, respectively: 

dP({J^})= d/i(J^), Jdn(Jij) = l, (7) 



p{ C )= n P ( ClJ ), ^ p(c iJ )=i. (8) 

(i,j),i<j <H,i=0A 



Depending on the temperature T, and on the param- 
eters of the probability distributions, dfi and p, the ran- 
dom model may stably stay either in the paramagnetic 
(P), in the ferromagnetic (F), or in the spin glass (SG) 
phase. In our approach for the F and SG phases there 
are two natural order parameters that will be indicated 
by m( F ) and m^ SG \ Similarly, for any correlation func- 
tion, quadratic or not, there are two natural quantities 
indicated by C' F ' and C^ SG \ and that in turn will be 
calculated in terms of and m^ SG \ respectively. To 
avoid confusion, it should be kept in mind that in our 
approach, for any observable O there are - in principle - 
always two solutions that we label as F and SG, but, for 
any temperature, only one of the two solutions is stable 
and useful in the thermodynamic limit. 

In the following, we will use the label o to specify that 
we are referring to the unperturbed model with Hamilto- 
nian Let mo(@Jo, Ph) be the stable magnetization of 
the unperturbed model with coupling Jo and in the pres- 
ence of a uniform external field h at inverse temperature 
(3. Then, the order parameters S=F,SG, satisfy 

the following self-consistent decoupled equations 

m< s > =mo( / 3J ( f ) , / 3J (s W E > + 0h), (12) 

where the effective couplings J^ F \ J^ SG \ Jq F ^ and Jq SG ' ) 
are given by 

f3J (F) =c J dn(J hj )t£inh(pj hj ), (13) 



The variables c^j € {0, 1} specify whether a "long-range" 
bond between the sites i and j is present (cjj = 1) or 
absent (cjj = 0), whereas the Jij's are the random vari- 
ables of the given bond For the Cjj's, we shall 
consider the following distribution 



(9) 



where c > 0. This choice leads in the thermodynamic 
limit N — > oo to a number of long range connections per 
site distributed according to a Poisson law with mean 
connectivity c. 

When we will need to be specific, for the Jij's we will 
assume either the distribution 



d£ ( Ji,j) 



$ (Ji,j - J) , 



(10) 



or 



dn{Ji,j) 
dJi j 



pd (Jij - J) dJu + (1 - p)S (Jij + J) , (11) 



to consider ferromagnetism or glassy phases, respectively. 
InEq. (TFrjpG [0,1]. 



dfi( Jij) tanh 2 (/3Jij), 



and 



r(F) _ 7 
J o — J o, 



P4 SG) = tanh _1 (tanh 2 (/3J )). 



(14) 



(15) 



(16) 



For the correlation functions , E=F,SG, we have 

= Q,(/?jf \ /?J< S W S ) + f3h) + (i) , (17) 

where Cb(/3Jo, (3h) is the correlation function of the 
unperturbed (non random) model. For the corrective 
0(l/iV) term in Eq. (jTTJ) we remind the reader to Eq. 
(31) of part I. Let us indicate by C' 1 ' and the aver- 
ages and the quadratic averages over the disorder of the 
correlation function of degree, say, k. and C^ 2 \ are 
related to C (F) and C (SG) , as follows 

C« = C (F) , inF, (18) 
C* (1) = 0, /codd, inSG, (19) 
C (1) = C (SG) , /ceven, in SG, (20) 
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and 



= (c( p )) 2 

C (2) = f C (SG)\'' 



in F. 



in SG. 



(21) 
(22) 



Among all the possible stable solutions of Eqs. (fl"2"l) . in 
the thermodynamic limit, for both £=F and S=SG, the 
true solution or leading solution, is the one that 

minimizes : 



= min L (s) (to) , (23) 
m£[— 1,1] 



where 



foiflJo, (3h) being the free energy density in the thermo- 
dynamic limit of the unperturbed model with coupling 
Jo and in the presence of an external field h, at inverse 
temperature (3. A necessary condition for a solution to/ s ' 
to be the leading solution is the stability condition: 



Xo jf ) , P J (E) m (S) + 0h) J (E) 



< 1, 



(25) 



where Xo(/?Jo, (3h) stands for the susceptibility xo of the 
unperturbed model divided by (3. 

For the localization and the reciprocal stability be- 
tween the F and SG phases we remind the reader to Sec. 
HID of part I. We recall however that, at least for lattices 
Co having only loops of even length, the stable P region 
is always that corresponding to a P-F phase diagram, so 
that in the P region the correlation functions must be 
calculated only through Eqs. (fl~8|) and pTjl . 

(X) 

If Jo > 0, the inverse critical temperature pc is solu- 
tion of the following exact equation 

Xo (/3f>jf\0)/3f)j( s ) = l, f <flW, (26) 

where f3 cQ is the inverse critical temperature of the un- 
perturbed model with coupling Jq . The constrain in 
Eq. (j26|) ensures the unicity of the solution. 

We end this section by stressing that this method is 
exact in all the P region and provides the exact critical 
surface, behavior and percolation threshold, and that, 
in the absence of frustration, the order parameters to/ s ' 
become exact also in the limit c — > + , in the case of 
second-order phase transitions, and in the limit c — > oo 
(see Sec. IIIC of part I). Note also that the order param- 
eters toA s ), and then the correlation functions, are by 
construction always exact in the zero temperature limit. 



IV. SMALL WORLD IN d = DIMENSION 

A. The Viana-Bray model 

As an immediate example, let us consider the Viana- 
Bray model. It can be seen as the simplest small-world 



model in which N spins with no short-range couplings 
(here Jo = 0) are randomly connected by long-range con- 
nections J (possibly also random). Note that formally 
here Co has dimension do = 0. Since Jo = 0, for the 
unperturbed model we have 



f3f (0,l3h)= log [2cosh(/%)] 



m (0,/3ft.) = tanh(/3/i), 



(27) 



(28) 



X~ (l3Jo,(3h) = l-tanh 2 (/3/i)|^ =0 = l, (29) 

It is interesting to check that the first and second deriva- 
tives of xo m h = are null and negative, respectively. 
In fact we have 

d 



and 



■^Xo(0,/%) = -2tanh 2 (/3/i) 
x [l-tanh 2 (/3/i)] |^ = o = 0, 



o2 

r X"o(0,j8ft) = -2 [l-tanh 2 (/%)] 2 



(30) 



4tanh 2 (^/i) [l - tanh 2 (^/i)] |^ =0 = -2. (31) 

Applying these results to Eqs. (|12j|16|) we get immedi- 
ately the self-consistent equations for the F and the SG 
magnetizations 



my' — tanh 
TO ( SG ) = tanh 



m (F) c / dA«tanh(/3J) 



m (SG) c / c^tanh 2 (/3J) 



(32) 



, (33) 



and the Viana-Bray critical surface 

(^tanh(/3( F) J) = 1, (34) 
c / ^tanh 2 (/3< SG) J) = 1. (35) 



On choosing for d[i a measure having average and vari- 
ance scaling as 0(l/c), for c cx V, we recover the equa- 
tions for the Sherrington-Kirkpatrick model already de- 
rived in this form in Q and Q. In these papers Eqs. 
(|32ti35[) were derived by mapping the Viana-Bray model 
and, similarly, the Sherrington-Kirkpatrick model to the 
non random fully connected Ising model. In this sense 
it should be also clear that, at least for /3 < f3 c and zero 
external field, in the thermodynamic limit, the connected 
correlation functions (of order k greater than 1) in the 
Sherrington-Kirkpatrick and in the Viana-Bray model 
are exactly zero. In fact, in the thermodynamic limit, 
the non random fully connected model can be exactly re- 
duced to a model of non interacting spins immersed in 
an effective medium so that among any two spins there 
is no correlation. Such a result is due to the fact that, in 
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these models, all the N spins interact through the same 
coupling J/N, no matter how far apart they are, and the 
net effect of this is that in the thermodynamic limit the 
system becomes equivalent to a collection of N non in- 
teracting spins seeing only an effective external field (the 
medium) like in Eq. (28) with (3h — ► /3Jmo. 

In the limit j3 -> oo, Eqs. (32) and (33) give the fol- 
lowing size (normalized to 1) of the giant connected com- 
ponent 



we get 



= tanh(m < - F - ) c), 
- tanh(m( SG )c). 



(36) 
(37) 



These equations are not exact, however they succeed 
in giving the exact percolation threshold c = 1. In fact, 
concerning the equation (3(5) f° r the F phase, the exact 
equation for m/ ¥ ^ is (see for example Q and references 
therein) 

l- m (F) = e - cm(F \ (38) 
which, in terms of the function tanh, becomes 

WF) + MF))2 tannic), (39) 



2-m( p ) + (m( p )) 2 

so that Eqs. (56") and (35) are equivalent at the order 
0(m^). We see also that, as stated in the previous 
Section, Eqs. (35) and (35) become equal in the limits 
c — > and c — > oo . 



B. Gas of Dimers 

Let us consider for £ a set of 2N spins coupled 
through a coupling Jo two by two. The expression "gas 
of dimers" stresses the fact that the dimers, i.e. the cou- 
ples of coupled spins, do not interact each other. As a 
consequence, the free energy, the magnetization, and the 
susceptibility of the unperturbed model can be immedi- 
ately calculated. We have 



1 



PfotfJo, Ph) = - log [2e? Jo cosh(2/%) + 2e 



(40) 



m (pj ,f3h) 



M s i n h(2/3ft,) 



e^ J ° cosh(2/3h) + e-f 3Ja ' 



(41) 



2e^°+2cosh(2/3fe) 

Xo{PJo,Ph) = — — — — 2\0h=o 

[e^ J ° cosh(2/3/i) + e~P Ja ] 



(42) 



cosh(/3J ) ' 

Let us calculate also the second derivative of xo ■ From 

■^Xo(pJo,ph)=4awh(ph) 
ph 

e -/3Jo _ 2e 3/3Jo _ e /3J COS h(2/3/l) 

[ e £ J ° cosh(2/?/i) + e"^] 3 



sinh(/3J ) 



a 3/3J 



[cosh(/3J )r 



(44) 



We note that, as expected, the second derivative of % in 
h = 0, for Jo > is always negative, whereas, for Jo < 
it becomes positive as soon as (3\Jo\ > log(v2). 

By using the above equations, from Sec. Ill we get 
immediately the following self-consistent equation for the 
magnetizations 



tanh(2/3j( s W s ) + 2/3h) 
1 + e -20J o sech(2/3 J( s )m( s ) + 2/3h) ' 



(45) 



and - at least for Jo > - the equation for the critical 
temperature 



cosh(M E) jf } ) 



-/3f)j( E ) 



(46) 



As it will be clear soon, this model lies between the 
Viana-Bray model and the more complex do = 1 dimen- 
sional chain small-world model, which will be analyzed 
in detail in the next section. Our major interest in this 
simpler gas of dimer small-world model is related to the 
fact that, in spite of its simplicity and do=0 dimension- 
ality - since the second derivative of xo may be positive 
when Jo is negative - according to the general result of 
Sec. IIIB of part I - it is already able to give rise to also 
multiple first- and second-order phase transitions. 



V. SMALL WORLD IN d =l DIMENSION 

In this section we will analyze the case in which £ is 
the g?o— 1-dimensional chain with periodic boundary con- 
ditions (p.b.c). The corresponding small-world model 
with Hamiltonian (2) in zero field has already been an- 
alyzed in Q by using the replica method. Here we will 
recover the results found in for f3 c and will provide 
the self-consistent equations for the magnetizations 
and m( SG ) whose solution, as expected, turns out to be in 
good agreement with the corresponding solutions found 
in Q for c small. It will be however rather evident how 
much the two methods differ in terms of simplicity and 
intuitive meaning. Furthermore, we will derive also an 
explicit expression for the two-points connected correla- 
tion function which, to the best of our knowledge, had not 
been published yet. Finally, we will analyze in the detail 
the completely novel scenario for the case Jo < which, 
as mentioned, produces multiple first- and second-order 
phase transitions. 

In order to apply the method of Sec. Ill we have to solve 
the one dimensional Ising model with p.b.c. immersed in 
an external field. The solution of this non random model 
is easy and well known (see for example l 9]). For the 
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free energy density, the magnetization and the two-points whereas for Jo < we have 
connected correlation function we have 

- (3f ((3J ,f3h) = log (Ax) , (47) dp^* o(/3J °' Ph) ~ for e ~ W ° > 3 ' (56) 



m, e pJo sinh(/3ft) 
m (/3 J , /3ft) = r , (48) 

[eWoginh 2 ^) + e- 2 ' 9J °] 2 



C ((3J ,(3h; \\i - j\\ Q ) = (cr t a j } - (a) 



X^\H-3\\ 
Ai, 



sin 2 (2^)(-^) ,(49) 



where the factor y> is defined by 



cot(2v3) = e 2/3J ° sinh^ft), 0< p < ~, (50) 

||i — j'Ho is the (euclidean) distance between i and j, and 
Ai and A2 are the two greatest eigenvalues of the matrix 
appearing in the transfer matrix method, whose ratio is 
given by (in the following expression the numerator and 
the denominator are A2 and Ai, respectively) 



A 2 eP Jo cosh(/3ft) - [e 2/3J ° sinh 2 (/3ft) + e - 2(3J °] 



Ai 



eP J ° cosh(/3ft) + [e 2 P J o sinh 2 (/3ft) + e" 2 ^] ' 



-.(51) 



Let us calculate \o and its first and second derivatives. 
From Eq. (|4*5|) we have 



Xo(/3Ja,fih) = 



-f}J cc . s h(/3ft) 



(52) 



We see therefore that, according to Sec. IIIB of part 
I, when Jo < for /3\J \ > log(3)/4 = 0.1193... we have 
a first-order phase transition. 

From Eqs. (fT2| and (|48|) . for the magnetizations 
and m^ SG ' at zero external field we have 



/?? 



s/ 3 J o sinh(/3J( E )m( s )) 



eW sinh 2 (/3J( s )m( s )) + e" 2 ^ 



(S) 2 



-• (57) 



From Eqs. (|26|) and (f52|) we see that a solution to( s ) 
becomes unstable at the inverse temperature (3 C given 
by 



e 2 ^' J o E Vf ' J (S) = 1. 



(58) 



For Jo > the above equation gives the exact P-F and 
P-SG critical temperatures in agreement with [8J. When 
Jo < - unless the transition be of second-order - Eq. 
(|58f for £ =F does not signal a phase transition. In 
general, as Jo < the P-F critical temperature must be 
determined by looking at all the stable solutions m of 
the self-consistent equation |57|) and by choosing the one 
minimizing the effective free energy L^(m) of Eq. 1(23)1 . 

Finally, for the two-point connected correlation func- 
tion, from Eqs. l)T?)). 1(45)1 and ((50]). we have 



d 



Xa(0J O> ph)=sinh(f3h) 



d(3h" 

[ e -3/3j _ 2e PJo cosh 2 ((3h) - e 13 ^] 



(53) 



C (S) (||i-j|| ) = sin 2 (2^ s ))e-H^llo/« <S) , (59) 



where 



2^ E )= cot" 1 e 2/3J o E) sinh(/3j( s W E )) , (60) 



2 XoC#Jo, Ph) = cosh(/3ft) 



d{(3h) 

[ e -3/3J _ 2eP J ° cosh 2 (/3h)-e 0Jo } t , 



r e 2/3Jo sin h 2 (/3/i) + e" 2 ' 3J «] 



0{(3h) 2 {$A) 



From Eq. ((54)) we see that for Jo > and any (3 we 
have, for sufficiently small ft, 



if 



d{(3h) 



-Xa{0Jo,Ph) < cosh(/3ft) 



[1 - 3eP J °] 



[e 2 P J °smh 2 (/3h) +e~ 2 P J °] 



- + 0(/3ft) 2 < 0,(55) 



and the correlation length is given by perform- 

ing the effective substitutions (3J — > /3Jq ' and /3ft — + 
/3j( s W s) in the logarithm of Eq. (J5T)) . 

Note that Co(/3Jo, /3ft) is even in /3ft, so that C(—m) = 
C(m). 

Near the critical temperature we have 



sin(2(^ 



e 2/3J^ E) m (S) 



+ 0(m^) 4 (61) 



and 
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log[tanh(/3 jf >)] - {fjJ m 

4sinh(/3J^ ) ) 



( e /"f > + ga/JjW) tanh(/3jf } ) + e 3 ^' - e^'] + 0(m( s >) 4 .(62) 
I 



One second-order phase transition 



No phase transition 




FIG. 1: Magnetization (solid line), and curves of stability 
(dashed and dot-dashed lines) for the case c = 0.5, Jo = 1 
and J = 3/5/. Here T c = 1.687. 
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0.2 - / 

" / 



5 6 7 



FIG. 3: Curves of stability (dashed and dot-dashed lines) for 
the case c = 5, Jo = —1 and J = 1. 



One second-order phase transition 




FIG. 2: Magnetization (solid line), and curves of stability 
(dashed and dot-dashed lines) for the case c = 10, Jo = 0.25 
and J = 1/c. Here T c = 1.419. 



No phase transition 



0.2 - ' 



4 5 6 7 



FIG. 4: Curves of stability (dashed and dot-dashed lines) for 
the case c = 5.828, Jo = —1.4 and J = 1. 



According to the general result (see Eqs. (51)-(55) of 
part I), we see that the correlation length remains finite 
at all temperatures. 

In Figs. 1-10 we plot the magnetization 
(solid line), xo (/? (F) ^ F) , o)P {F) J (F) (dashed line), and 

Xo (/? (F) 4 F) :/ ?j(F)m(F) )/ 3(F)j(F) (dot-dashed line) for 
several cases obtained by solving Eq. (|5T[) numerically 
with E=F, and by choosing the stable solution minimiz- 
ing l/ F )(m) (see Eq. (|24p). In all these examples we 
have chosen the measure (|10l) . Figs. 1 and 2 concern 
two cases with Jo > so that one and only one second- 
order phase transition is present. The input data of these 
two cases are the same as those analyzed numerically in 
Q (note that in the model considered in [8], the long 



range coupling J is divided by c). As already stated in 
Sec. Ill, the self-consistent equations become exact in 
the limit c — ► 0, for second order phase transitions, and 
in the limit c — > oo. Therefore, for the magnetization, by 
comparison with Q, in Fig. 1 and 2, where c is relatively 
small and big, respectively, we see good agreement also 
below the critical temperature. 

Figs. 3-10 concern eight cases with Jo < 0. As ex- 
plained above, the critical behavior and the localization 
of the critical temperatures is more complicated when 
Jo < 0. In particular, given | Jo|, if c is not sufficiently 
high the solution = remains stable at all tem- 

peratures and if it is also a leading solution, no phase 
transition occurs. Let us consider Eq. For J < 

the lhs of this equation has some maximum at a finite 
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1^=3.641, two first-order phase transitions 



1^=1.820, first- and second-order phase transition 
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FIG. 5: Magnetization (solid line), and curves of stability 
(dashed and dot-dashed lines) for the case c = 5.5, Jo = — 1, 
and J = 1. Here T c i = 1.02 and T c2 = 2.27. 

1^=3.277, two first-order phase transitions 



FIG. 6: Magnetization (solid line), and curves of stability 
(dashed and dot-dashed lines) for the case c = 5.5, Jo = —0.9 
and J = 1. Here T cl = 0.85 and T c2 = 2.81. 



value f3 given by 
0J 



log 



1 + 5{r) 



1 - 5{r) 



where r = I Jq\/ J, and we have introduced 



S(r) = \J\ + r 2 - 



(63) 



(64) 



Hence, we see that a sufficient condition for the solution 
ttx ' = to become unstable is that be 



1 + 5{r) 
1 - S{r) 



8(r) > 1. 



(65) 



Note that the above represents only a condition for the 
instability of the solution — 0, but the true solution 
is the one that is both stable and leading. In fact, when 
Jo < 0, a phase transition in general may be present also 
when Eq. (|6"5|) is not satisfied and, correspondingly the 
possible critical temperatures will be not determined by 

Eq. (EHD- 

In Fig 3 we report a case with J = 1, Jq = — 1 (r = 1) 
and a relatively low value of c, c = 5, so that no phase 




FIG. 7: Magnetization (solid line), and curves of stability 
(dashed and dot-dashed lines) for the case c = 6, Jo = —0.5 
and J = 1. Here T c i = 0.35 and T C 2 = 4.87. T c i corresponds 
to a second-order phase transition. 

T,=0.728, first- and second- order phase transitions 




FIG. 8: Magnetization (solid line), and curves of stability 
(dashed and dot-dashed lines) for the case c = 4, Jo = —0.2 
and J — 2. Here T a \ = 0.22 and T C 2 = 7.55. T C 2 corresponds 
to a second-order phase transition. 



transition is present. Similarly, in Fig. 4 we report again 
a case in which no phase transition is present due to the 
fact that here r is relatively big, r = 1.1. It is interesting 
to observe that for r = 1 Eq. (|6"5|) requires a value of 
c greater than the limit value c = 3 + 2v2 = 5.8284.... 
In both of Figs. 5 and 6 we report a case in which Eq. 
(IBT))) is still not satisfied, but nevertheless two first-order 
phase transitions are present. In both of Figs. 7 and 
8 we have one first- and one second-order phase transi- 
tion. In both of Figs. 9 and 10 we have two second-order 
phase transitions. The critical temperature of second- 
order transition can be determined also by Eq. (f58|) . In 
the top of Figs. 5-10 we write the discriminant temper- 
ature T» = 4|Jo|/log(3) below which a phase transition 
(if any) is first-order (see Eqs. (|5"4")) - (l5"6"l) ) (see Sec. IIIB 
of part I). 

Finally in Figs. 11 and 12 we plot the 
spin glass order parameter m' SG ' (solid line), 

Xo (/? (SG) Jo SG) ,0)/3( SG )j( SG ) (dashed line), and 
Xo (/3( SG )J^ SG) ,/?J( SG )m( SG ))/?( SG )j( SG ) (dot-dashed 
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T s =2.185, two second-order phase transitions 



0.2 - / 




FIG. 9: Magnetization (solid line), and curves of stability 
(dashed and dot-dashed lines) for the case c = 1.6, Jo = —0.6 
and J = 7. Here T cl = 2.58 and T c2 = 7.55. 

T # =1.820, two second-order phase transitions 



/ 




FIG. 10: Magnetization (solid line), and curves of stability 
(dashed and dot-dashed lines) for the case c = 1.4, Jo = —0.5 
and J = 10. Here T cl = 3.00 and T c2 = 9.34. 

line) obtained by solving Eq. (|57p numerically with 
£=SG. In these two examples we have chosen the mea- 
sure (fTTj) and, for c, J, and Jo, we have considered the 
same parameters of Figs. 1 and 2 of the ferromagnetic 
case. 

Note that, unlike the P-F critical surface, the P-SG 
critical surface does not depend on the parameter p enter- 
ing in Eq. (fTTj) . For the reciprocal stability between the 
P-F and the P-SG critical surface we remind the reader 
to the general rules of Sec. HID of part I (see cases (1) 
and (3)) which, for J > 0, reduce to the results reported 
in Sec. 6.1 of the Ref. @]. Here we stress just that, if 
Jo > 0, for p < 0.5, only the P-SG transition is possi- 
ble. However, when Jo < and c is not sufficiently large, 
the SG phase may be the only stable phase even when 
p = 1. In fact, although when Jo < the solution 
may have two P-F critical temperatures, in general, if the 
P-SG temperature is between these, we cannot exclude 
that the solution m^ SG ^ starts to be the leading solution 
at sufficiently low temperatures. In Figs. [T3landfT4l on 
the plane (T, c) , we plot the phase diagrams correspond- 
ing to the cases of Figs. and QUI respectively. Theses 
phase diagrams are obtained by solving Eq. (|26p suppos- 



Spin glass second-order phase transition 




0.5 1 1.5 2 

T 

FIG. 11: Spin glass order parameter (solid line), and curves 
of stability (dashed and dot-dashed lines) for the case c = 0.5, 
Jo = 1 and J = 3/5/c. Here T c = 1.130. 



Spin glass second-order phase transition 




0.5 1 1.5 2 

T 



FIG. 12: Spin glass order parameter (solid line), and curves 
of stability (dashed and dot-dashed lines) for the case c = 10, 
Jo = 0.25 and J = 1/c. Here T c = 0.424. 



ing that here, as in the cases of Figs. [9] and [101 where 
c = 1.4 and c = 0.5, respectively, the P-F transition is al- 
ways second-order. We plan to investigate in more detail 
the phase diagram in future works. 



VI. SMALL- WORLD SPHERICAL MODEL IN 
ARBITRARY DIMENSION do 

In this section we will analyze the case in which the un- 
perturbed model is the spherical model built up over a Jo- 
dimensional lattice Co (see Q and references therein) [J] . 
In this case the cr's are continuous "spin" variables rang- 
ing in the interval (— oo,oo) subjected to the sole con- 
strain YV erf = N, however our theorems and formal- 
ism can be applied as well and give results that, within 
the same limitations prescribed in Sec. Ill, are exact. 

Following [9(, for the unperturbed model we have 

- /3/o(/3J , (3h) = I log + <$> (13 Jo, ph, z) , (66) 
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Phase diagram 




FIG. 13: Phase diagram for the case considered in Fig. |9]with 
the measure of Eq. (|TD|) . 

Phase diagram 




FIG. 14: Phase diagram for the case considered in Fig. [10] 
with the measure of Eq. (I10|l . 



m {PJ ,/3h) 



0h 
2/?J z' 



where 



1 ( 3h) 
(j) ((3J ,f3h, z) = [3J d + (3J z - -g(z) + i^-L 

2 4/jJqz 



(67) 



, (68) 



1 



2tt 



2n 



W l ° JO JO 
x log [d + z — cos(wi) 



dui . . . duj do 

. . - cos(u do )] , (69) 



and z = z((3Jq, f3h) is the (unique) solution of the equa- 
tion d z cj> (PJo, flh, z) = 0: 



PJo 



(my 



i 



4/3 Jo z 2 T 
from which it follows the equation for tuq 



(3J (1 - ml) = -g 



1 , ( @h 



2 a VWoTOo, 
The derivative g' can in turn be expressed as 

/>OC 

g'(z) = / e -*(*+*) [Jo (»t)]* dfc, 
Jo 



(70) 



(71) 



(72) 



jTo (it) being the usual Bessel function whose behavior 
for large t is given by 



Jo (it) 



1 + 



(73) 



The critical behavior of the unperturbed system de- 
pends on the values of g'(z) and g"(z) near z = 0. It 
turns out that for do < 2 one has g'(0) = oo and there 
is no spontaneous magnetization, whereas for do > 2 one 
has <?'(0) < oo and at h = the unperturbed system 
undergoes a second-order phase transition with magneti- 
zation given by Eq. (|7Tj) which, for j3 above f3 c o, becomes 



(74) 



m (/3Jo,0) = 4/1 -^r 




where the inverse critical temperature /3 c q is given by 



Aa Jo = \g' (0) • 



(75) 



Furthermore, it turns out that for do < 4 one has 
g"(0) — oo, whereas for do > 4 one has g" (0) < oo. 
This reflects on the critical exponents a, 7 and 5, which 
take the classical mean-field values only for do > 4. 

According to sec. Ill, to solve the random model - 
for simplicity - at zero external field, we have to per- 
form the effective substitutions p Jo — > /3Jq and (3h — » 
/3J( E W S ) 



in the above equations. From Eqs. (|67j) , 1701 



and (|7ip . we get immediately: 



/3J( S ) 



2/3J 



(£)> 



the equations for inverse critical temperature (3, 



(E) T(S) _1 , ( fr>J™ 
" 2 5 I 2/3 f)jf ) 



P^J'o 

and the magnetizations mP" 



(76) 



(77) 



777. 




2 /"o 

£</3, 



/3J( S > 



> 



1:1 



/8>/9, 



(78) 



Note that, as it must be from the general result of Sec. 
IIIB of part I, unlike the unperturbed model, as soon 
as the connectivity c is not zero, Eq. (|TT[) has always a 
finite solution pc , independently on the dimension do- 
In fact, one has a finite temperature second-order phase 
transition even for do — > + where from Eq. (|T2")) we have 



1 



9'(z) = - 



(do = 0) 



(79) 



so that the equations for the critical temperature ([77 
become 



/5f)j< 5 



1, 



(do = 0) 



(80) 
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which, as expected, coincide with Eqs. (|34|) and (|35[) of 
the Viana-Bray model. 

Similarly, unlike the unperturbed model, in the ran- 
dom model all the critical exponents take the classical 
mean- field values, independently on the dimension do- 
In the specific case of the spherical model, this behavior 
is due to the fact that g'(z) and g"(z) can be singular 
only at z = but, as soon as the connectivity c is not 
zero, there is an effective external field (3J^m^ so that 
is not zero. For the critical behavior, the dependence 
on the dimension rf reflects only in the coefficients, not 



on the critical exponents. In particular, concerning the 

argument of the square root of the rhs of Eq. l[78]). by 

expanding in the reduced temperature , for | <C 1 
we have 

1 i_ , f ^J^ ) 

£>(£)£(£) +0(< (£))2 ; (81) 

where 



.(2/3f ) 4 F) ) 2 



dtt(Ju) (l - tanh 2 03f > Jy)) /9f } - c ^ ^(J^) tanh(/3< p > , (82) 



1 



(SG) 7 (SG) 



_ 4 tanhjM^Jo)^Jo + „ 
l + tanh 2 (^ SG) J ) 



(2^ SG) jf G) )2 



2c 



d/i(Jij) (l - tanh 2 (4 SG ) Jij)) tanhCSf 30 ) J id )0i SG > J;, 



-4c / dtffj- )tanh 2 (/3( SG )j- ■) tanh(^ 5G) J )/# G) J 

J M M) ^ (l + tanh 2 (/# G > J o) ) M SG ) J(SG) 



(83) 



so that from Eqs. (|78[) and (fSTj) for the critical behavior 
of the magnetizations we get explicitly the mean-field 
behavior: 

\ 0, ^ > 0. K ' 

VII. CONCLUSIONS 

We have applied a novel method to solve analytically 
several small- world models of interest defined over an un- 
derlying lattice Cq of dimension do = 0, l,and oo, corre- 
sponding to an ensemble of non interacting units (spins, 
dimers, etc.), the one dimensional chain, and the spher- 
ical model, respectively. The long-range couplings may 
also have an additional own disorder leading, in partic- 
ular, to spin glass phases. The simplicity of the method 
allows us to find very easily (just a few line equations) 
the critical surfaces and the order parameters. In our 
method the former are exact whereas the latter provide 
an approximation to the exact order parameters which, 
in the absence of frustration, become exact in the limit 
c — ► 0, when the transition is of second order, and in the 
limit c — > oo. In particular, we have studied in detail the 



small-world defined over the one dimensional chain with 
positive and negative short-range couplings, showing ex- 
plicitly how, in the second case, multicritical points with 
first- and second-order phase transition arise. Finally, the 
small- world spherical model - an exact solvable model (in 
our approach) with continuous spin variables - has pro- 
vided us an interesting testing bench case to study what 
happens as do changes continuously from to oo. As 
expected from general grounds, unlike the non random 
version of the model, the small- world model presents al- 
ways a finite temperature phase transition, even in the 
limit do — + + . This latter result - besides to be consis- 
tent with what we have found in the do = 0-dimensional 
discrete models - is easily and physically explained by 
the method. In fact, it consists on mapping the small- 
world model (a random model) to a corresponding non 
random model (no long-range bonds) but immersed in 
an effective uniform external field which is active as soon 
as the added random connectivity c is not zero (see Eqs. 

(3D-CE1 ). 

In this paper we have considered small-world models 
which - by the method - turn out to be the simplest 
ones (cases for which the unperturbed model is analyti- 
cally feasible) . Many interesting variants of these models 
can be as well considered and are still analytically solv- 
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able by the method. However, the method can be also 
applied numerically to study other more complex small- 
world models for which the corresponding unperturbed 
model is not analytically solvable. In fact, the numerical 
complexity in solving such a small-world model, will be 
still as feasible as a non random model immersed in a 
uniform external field. 
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